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tivity is calculated for a superconductor. Sharp structure, characteristic of 
superconductivity with an order parameter with s-wave symmetry, appears in 



^ , the BCS limit as a minimum at a frequency equal to twice the gap value. This 

^ structure scales with temperature but gets progressively smeared and shifted 

^ , as impurity scattering is increased. The relationship between low frequency 

H : 

results and the zero frequency limit is investigated. Experimental results on 
Bai^xKxBiO^ are also presented. The low frequency imaginary part of the 
conductivity displays a minimum at 2A ~ 12 meV, and provides unequivo- 
cal evidence of an s-wave superconducting order parameter. Strong coupling 
(Eliashberg) results show similar trends. Using this formulation we find that 
the electron-phonon coupling in Bai^xKxBiO^ must necessarily be small, 
with coupling constant A ~ 0.2, in agreement with conclusions drawn from 
measurements of the real part of the conductivity. Thus Bai-xKxBiO^ is an 
s-wave superconductor that is not driven by the electron-phonon interaction. 
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I. INTRODUCTION 

Optical conductivity data in the infrared can give important information on the proper- 
ties of the low energy excitations present, and, in particular, on the scattering rate [|I|. In 
Drude theory, the width of the real part of the conductivity gives the inverse of the scat- 
tering time, which can be a combination of elastic impurity scattering and inelastic (e.g. 
phonon-mediated) scattering. In the latter case, a new process, phonon-assisted absorption, 
becomes possible P] . This process does not exist with impurity scattering alone, and makes 
an additional contribution to both the real and imaginary parts of the conductivity. 

In the past most discussions of the infrared conductivity in an electron-phonon system 
have concentrated on the real part of the conductivity [|l|. Very few results are available 
on the imaginary part which nonetheless follows uniquely from the real part through a 
Kramers-Kronig relation. Recently experimentalists have effectively exploited the relation- 
ship between the inverse square of the penetration depth and the zero frequency limit of the 
imaginary part of the conductivity {limy^Qva2{^) = |^a2^5 where we use the convention 
that a = o"i — ia2i and A here is the penetration depth), to obtain information on this 
important quantity from the far-infrared optical conductivity . 

The purpose of this paper is to study the frequency dependence of ua2{i')- We start 
with the BCS model and then generalize to an Eliashberg strong coupling system. The 
aim is to understand how the more complicated effects due to the interactions with phonons 
modify the BCS results. The necessary theory is summarized in section II. Results of BCS 
theory are presented in section III. It is found that in an s-wave superconductor there will 
be a distinct minimum in z/(T2(z/) at twice the gap value. With the addition of impurities the 
minimum broadens and shifts slightly. The temperature dependence of the minimum tracks 
the gap temperature dependence. We also study how finite frequency results are related 
to the zero frequency limit. Results are presented in section IV for the strong coupling 
case. The most significant difference between the Eliashberg and BCS approaches occurs 
in the pure case. In this limit, BCS theory no longer contains any scattering mechanisms, 
whereas inelastic scattering remains in an Eliashberg formalism, and we can study the effect 
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of this additional source of scattering. In section V, the specific case of Bai-^K^BiOs (with 
X = 0.4) will be considered. Comparison with infrared data will lead us to conclude that 
this system contains an isotropic s-wave order parameter, but is not an electron-phonon 
superconductor. The inelastic scattering observed in the conductivity at room temperature 
is inconsistent with a mass enhancement parameter of order unity, which is required for an 
electron-phonon superconductor with Tc 30 K. A summary is provided in section VI. 

II. FORMALISM 

Expressions for the optical conductivity of an electron-phonon superconductor, excluding 
vertex corrections for the electron-phonon interactions, are now well established |p|^!^ and 



require solutions of the Eliashberg equations [H for the pairing function (f){uj) and renor- 
malized frequency uj{uj). For an isotropic s-wave electron-phonon superconductor, these 
equations take the form |^ 



m=—oo 

+171 r dv a^F{v) I [N{v) + f{v - uo)] 



m=-oo MUj'^{iuJra) + 4>m 



^uj'^{lj — u) — (p{ljJ — u) 



and 



+ m + /(. + .)]^^h±^^_| (1) 



(jj[UJ) = 1 + tTTl A[UJ — lUJm) = 

+Z7r dv a^F{v) \ [N{v) + f{v - uo)] 



U0[LO — V] 



^Lj'^{uJ — u) — (j)^{uj — u) 

+ [iVM + + }. (2) 

In these equations Nlu) and /(z/) are the Bose and Fermi distribution functions, respectively, 
and X{z) is the Hilbert transform of the electron-phonon spectral function a'^F{i'). fi*{ujc) 
is the Coulomb repulsion parameter with cutoff Uf.. A negative value for this parameter can 
be used to model some BCS attraction of unspecified origin. The required functions are first 



obtained on the imaginary axis at the Matsubara frequencies, i.e. uo = iujn = iTiT{2n — 1), 
with 0m = 4'{i^m) by setting the complex variable u in these equations to the Matsubara 
frequencies [|I6],|17|. Then the equations are iterated as written, with u set to a frequency on 



the real axis. 

From the numerical solutions of Eqs. (|T]-0), the conductivity follows from the paramag- 
netic response kernel, n(z/ + i6), given by i5|j6|J^^13| 

n(z/ + z5) = ■( — 1 + / (iwtanhf — \\hx{io,LO ^ v\ — h^iyi.io ^ v)] 

m [ JO 2 ^ ' 

+ J"" dujt&nhi ^^'^^ ''^ )[hl{uj,uj + v) + h2{u:,uj + v))^ (3) 

with the conductivity 

a{v) = -(u{v + i5) + —), (4) 

where e is the magnitude of the charge on the electron, m its mass and n the number of 
electrons per unit volume. The functions hi in Eq. (j^) are given by 

1 - N{UJI)N{U2) - P{UJI)P{U2) 



hi{uJi,UJ2 
h2{uJi,UJ2 

N{uj 
Piuj 



2{e{uJi) + e{uj2)) 
1 + N*{U^)N{U2) + P*{ui)P{uJ2) 

2{e{u2) - e*(a;i)) 
Cj{uj + i5) 
e{uj + i5) 



e{uLi + iS) 

e{u) = ^Cj'^{u} + i5) - (t)'^{u + i5). (5) 

The equations so far have been written for a clean system. To include normal impurity 
scattering in the Born approximation, an average scattering rate, 1/r, is used, and gives rise 
to an additional contribution to the pairing and renormalization functions: 

i:;{uj) ^ uju) + - , ^^^^ (7) 



Eq. (0) remains the same with impurity scattering. The modifications are all implicitly 
contained in the pairing and renormalization functions. Note that the gap parameter, 
A(u;) = (p{uj) I Z(u), remains the same, independent of the impurity scattering rate. 

III. THE BCS LIMIT 

The equations for the conductivity and pairing function discussed so far reduce in the 
BCS limit, thus simplifying the numerical work |T^^0|. The result is a universal function for 



the conductivity, dependent only on the impurity scattering rate. In the BCS limit uj{uj) is 
given by the bare frequency u (plus the impurity scattering term given in Eq. (0)), and the 
pairing function (p^uj) reduces to the BCS gap parameter, A (plus Eq. (P)). The temperature 
dependence of A(T) is given by solution of the BCS equation. As already mentioned, the 
clean limit is pathological within BCS theory. To allow some scattering mechanism to remain 
some nonzero impurity scattering is required. Note that in the remainder of this paper we 
present results for the conductivity normalized to ne'^/m, so the only parameter in the BCS 
limit is the impurity scattering rate, 1/r. 

Results are shown in Fig. 1 for h'a2{i') as a function of frequency for a BCS model with 
T, = 29 K (Bai^^K^BiOs). We use (a) 1/r = 1, (b) 1/r = 5, and (c) 1/r = 25 meV, with 
differently styled curves for various temperatures, as explained in the figure caption. The 
largest scattering rate we have used is representative of the Bai-xK^BiO^ compounds, as 
we will see later. 

For small scattering rate (Fig. (la)) the quantity i>(T2{i') displays a cusp-like minimum 
at a frequency equal to two times the temperature dependent gap parameter. As the tem- 
perature increases towards the size of the minimum decreases, and of course its location 
tracks the energy gap. For 1/r = 1 meV the cusp-like minimum occurs at perhaps too low 
a frequency for temperatures near to be observed with optical infrared techniques, but at 
least in principle this minimum can be used to measure the gap temperature dependence. 
Note that the cusp in the imaginary part of the conductivity simply reflects the sharp onset 
in the real part of the conductivity which occurs at 2A(T). For larger impurity scattering 



5 



rates the minimum broadens somewhat. 

Another important feature of our results is the hmit as the frequency goes to zero, 
which as stated earher, gives the London penetration depth. We have verified that our low 
frequency conductivity results extrapolate to the correct penetration depth, as computed 
from an imaginary axis formulation P,pT|: 

As is apparent from the figures, particularly for low scattering rate, this limit can only 
be probed with very low photon energy, certainly too low for the usual infrared range in 
conventional electron-phonon superconductors. 

As the impurity scattering rate increases the low frequency conductivity i^o"2(z/) deviates 
significantly from unity. This is expected since the penetration depth increases with impurity 
scattering in a BCS model. The effect is very pronounced in Fig. (Ic) where it is also seen 
that the minimum due to the gap is considerably broadened due to the impurity scattering. 
To see the effect of impurity scattering on the low frequency conductivity more clearly we 
show in Fig. 2 results for various impurity scattering rates. All curves are calculated in the 
BCS limit and for reduced temperature, T/T^ = 0.1, which is essentially zero temperature. 
As the scattering rate increases, the cusp-like minimum broadens into a shallow minimum. 
For very high scattering rates, it is evident that extracting a gap value would be very difficult. 
The actual minimum is displaced to higher frequencies, and is no longer located at twice 
the gap. 

The results of Figs. (1-2) immediately suggest the following question: To what extent 
can finite frequency conductivity data be used to extract information about the zero fre- 
quency penetration depth ? We first address the question of temperature dependence. The 
dimensionless quantity z/cr2(z/) is plotted in Fig. 3 as a function of reduced temperature, for 
various frequencies, as given in the figure caption. Two impurity scattering rates are used, 
(a) 1/r = 1 meV to represent a fairly clean case, and (b) 1/r = 25 meV, to represent both 
a dirtier case and Bai^xKxBiO^. In Fig. (3a) the lowest frequency curve plotted {v = 0.2 
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meV) does indeed follow very closely the temperature dependence of the zero frequency 
penetration depth. Nonetheless, already at z/ = 1.0 meV (dotted curve) the temperature 
dependence of z/cr2(z^) deviates very markedly from that of the solid curve, and hence cannot 
be used to extract the temperature dependence of the penetration depth. As the impurity 
scattering increases the finite frequency results improve as is seen in Fig. (3b). Clearly the 
finite frequency results give a good qualitative estimate of the temperature variation of the 
penetration depth up to u = 2.5 meV. For frequencies of order twice the gap (e.g. u = 10 
meV), however, all correspondence with the penetration depth is lost. 

Can finite frequency data be used to extract the impurity dependence of the penetration 
depth ? In the BCS limit, with a = [||], 



1 1 \ 71 1 _i , , ^, 1 

sin (vl - a )( a <1 



A2 (T = 0) A2 (T = 0) j 2« aVT^ 



1 



f TT 1 , /a + Vtt^ — 1\ 1 , , 

— 

a;^ — 1 — Va'' — 1 ^ 



In the weak scattering limit this expression reduces to the more familiar form, 
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(10) 



A2(0) A2(0)l + f« 

This expression can be written in terms of the zero temperature coherence length, ^q? and 
the mean free path, using A = and vp = ^/t, where vp is the Fermi velocity: 

111 

Ap) TTff' ^^^^ 

In Fig. 4 we plot z/cr2(z/) at various frequencies as a function of the scattering rate, 1/r. 
Clearly, up to about 5 meV the deviations from the zero frequency curve, Eq. (P), are small. 
Significant deviations begin to occur at higher frequencies {u = 10 meV). Thus, one must 
be careful when using finite frequency data to infer information about the zero frequency 
penetration depth. 

IV. THE STRONG COUPLING CASE 

The results presented so far are based on BCS theory, where the interaction is modelled 
by a constant, structureless in momentum and frequency (except for a cutoff). One of the 



merits of such an approach is that the results are universal, but the disadvantage is that 
inelastic scattering is not included and so strong coupling effects in the superconducting state 
are not accounted for. Strong coupling results first require a full numerical solution to the 
Eliashberg equations (0,0). Two parameters are needed to characterize a particular material, 
the electron-phonon spectral function a'^F{uj), and the Coulomb repulsion, /i*. The latter 
is often fitted to either Tc or the gap, while the former is extracted through an inversion 
procedure from tunnelling data Often a'^F{uj) has been found to be given roughly by 
a constant times the phonon frequency distribution, F{uj). Thus, for most purposes, the 
difference between a'^F{uj) and -^(ci;) can be ignored and a single multiplicative parameter 
can be used to obtain a'^F{uj) from F{uj) when tunnelling data is not available. To be specific 
here, we will treat only the case of Bai^^KxEiO^. Its frequency distribution is known 



231 and is similar to calculated [E5|l values of a'^F(uj) as well as values determined from 



tunnelling p6[. Based on these three spectra then if Bai-xK^BiO^ is an electron-phonon 
superconductor, the value of the electron mass renormalization A = 2/ duja'^F{(jj) / uj needed 
to produce superconductivity with Tc ~ 30 K is A ~ 1. Because of difficulties with the 
inversion of the tunnelling data PB|, we will use a spectrum scaled from neutron scattering 



results with A = 1 as a representative intermediate coupling Eliashberg superconductor. In 
fact Sharifi et al. have argued, based on tunnelling measurements, that Bai-xK^BiO^ is not 
an electron-phonon based superconductor In the next section we will also argue that 
Bai^xK^BiO^ cannot be an electron-phonon superconductor and that A is more likely to be 
of order 0.2, a value much smaller than considered in the work previously quoted [p^-p6[|. 

In Fig. (5a) we show results for the product of the frequency times the imaginary part 
of the conductivity, i.e. z/cr2(z/), as a function of z/ for a strong coupling superconductor 
with A = 1 and = 29 K as discussed above. The curves are plotted in the clean limit, 
i.e. 1/r = 0. Scattering is still possible through inelastic scattering via the phonons, which 
are present at finite temperature. At low temperatures a very sharp Drude-like peak will 
be present in the normal state conductivity. This manifests itself in the imaginary part of 



the conductivity as a sharp inverted Drude-like minimum near the origin [^, as is visible 



in the lower left portion of Fig. (5a). The solid curve is for T/Tc = 0.99 while the dotted 
curve, which is almost indistinguishable from the solid curve, is for the normal state at the 
same temperature. This curve extends down to zero at the origin (not shown) and has a 
width of order 1 meV. Also shown are results at T/Tc = 0.5 (dashed curve) and T/Tc = 0.2 
(dot-dashed curve). These curves are very different from their counterparts in the previous 
section. The first thing to note is that sharp structure at z/ = 2A ^ 10.4 meV is completely 
absent. In this case there is no BCS reference, as in this limit the corresponding result would 
be a constant at value unity (see Fig. 1). Recall that in the real part of the conductivity 
absorption begins beyond 2A not in an abrupt fashion as in the Mattis-Bardeen limit of 
the theory, but smoothly, roughly following the phonon density of states since the phonons 
conserve the energy and account for the momentum of the created electron-hole pair |]8|,p9|. 
Since the absorption edge is not abrupt in cri(z/), we cannot expect the corresponding sharp 
structure in the imaginary part h'a2{i'); this is born out by our calculations shown in Fig. 
(5a). 

Another difference which is obvious upon inspection of Fig. (5a) is that structure is 
present in the imaginary part of the conductivity over a frequency range representative of 
the phonon energies in the a'^F{uj) spectrum. They occur here in much the same way as in 
the real part of the conductivity |^,[13|, and leave open the possibility for the determination 



of a^F^Lj) by infrared spectroscopy. 

Finally another clear difference occurs in the zero frequency limit, where, at low tem- 
peratures, the intercept is approximately 0.5, i.e. half of the London limit. This is due to 
phonon renormalization effects, which, for the most part amount to a 1 + A enhancement of 



the mass, as is obvious from Eq. (||) pO |. 



For non-zero impurity scattering two important things happen, as Fig. (5b) shows. A 
cusp-like minimum, conspicuous by its absence in the clean limit, is now present even with 
a small impurity scattering rate, 1/r = 2 meV. Secondly the penetration depth, given by 
the inverse of the zero frequency limit of the quantity plotted (see remark in the introduc- 
tion), has increased with the addition of impurities, as discussed in the previous section in 
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the BCS limit. The increase here is non-universal, i.e. it depends on the strength of the 
interaction, and somewhat less sensitively on the form of the underlying electron-phonon 
spectral function. 

To see the effect of coupling strength on the low frequency imaginary part of the con- 
ductivity, the quantity h'(T2{i') is plotted in Fig. 6 vs. frequency with (a) 1/r = 2, and (b) 
1/r = 25 meV. The curves plotted are for various coupling strengths as indicated, with the 
Coulomb repulsion /i* adjusted in each case so that Tc remains 29 K. We have used a small 
Coulomb repulsion to model Bai-xK^BiO^ with A = 1 so that for small A we use negative 
/i*'s, thus representing some additional attractive mechanism by which the electrons bind. 
In Fig. (6a) there is a clear lowering of the imaginary part of the conductivity with increas- 
ing coupling strength, while in Fig. (6b) the value of va2{i') is drastically reduced even in 
the BCS limit (A = 0). In both cases the depth and sharpness of the minimum associated 
with 2A decreases as the coupling strength increases. The minimum also shifts slightly to 
higher frequency, as expected since the zero temperature gap increases with increasing cou- 
pling strength. At higher frequencies phonon structure becomes more noticeable as coupling 
strength increases. 

In Fig. 7 we show similar results to those given in Fig. 5 but now with an impurity 
scattering rate equal to 50 meV. Upon comparison with Fig. 5, it is clear that the entire 
curve for va2{i') at low temperatures {T/Tc = 0.2) has been reduced systematically by the 
impurity scattering. Furthermore a gap minimum is clearly present, as Fig. (5b) would 
suggest. Phonon structure is still seen at higher energies. Also, the value of the penetration 
depth has been reduced significantly in a fashion that roughly obeys Eq. (|^). There are 
strong coupling corrections ||2l[ so that exact agreement is not expected. 



The higher temperature curves are also quite interesting and deserve some attention. 
The solid curve is for the superconducting case with T/Tc = 0.99 and falls almost exactly 
on the dotted curve which is the normal state result at the same temperature (differences 
are only apparent at very low frequency). In the normal state, the low frequency part of the 
real and imaginary part of the conductivity can be fit to a Drude-like contribution, as first 
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attempted by Dolgov et al. ||3T| and taken up more microscopically and rigorously recently 
by us |13[ . If both real and imaginary parts are fitted in a low frequency expansion of the 
normal state conductivity, we arrive at the result 

_ ne^ If 

O'Drude = :; ~Z r-T, (12) 



with 



m 1 -)- ;\ 1 — iuT 



and 



with 



f = ^ ^, (13) 

(r*(^)/(l + A*(^))) 



l + X{u) = ^^ (15) 



and 



l/r*{u) ^ V^ + ^Im^H . (16) 

Note that {Q{u)) = ^^sech '^{^)Q{uj) and explicit closed expressions for both these 
last quantities are given in our previous paper |[T3|| . 

In Fig. 8 we compare our Drude fits (dotted curves) with the full electron-phonon results 
(solid curves) in the clean limit as a function of frequency for various temperatures in the 
normal state. At low temperature, i.e. T/Tc = 1, the Drude fit at low frequency is very 
narrow, reflecting the fact that the electron-phonon scattering time becomes very long at 
low temperature. As the frequency increases the two curves begin to deviate very markedly. 
In units of ne^/m, z/cr2(z^) as modelled by the Drude form Eq. (0) will saturate at high 
frequencies to a value of 1/(1 + A), which is 0.43 for the chosen parameters. This saturation 
is seen very clearly in the plots. As the temperature is increased to T/T^ = 5, the inelastic 
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scattering increases significantly. The solid curve reflects this fact in that the low frequency 
variation is on a much larger scale, of order 60 meV. Also, the Drude-like curve (dotted) 
fits the full numerical solution (solid curve) over a larger frequency range and saturates at a 
larger value of 1/(1 + A), which has increased considerably with increasing temperature as 
described in our previous work Note that lim^^^oo z/(T2(z^) = 1 in the full calculation (in 
units of ne^/m). Thus, one could use i/a"2(z^) data to determine A and f independently from 
the low frequency fit and from the infinite frequency limit. Finally, the last set of curves is 
for T/Tc = 10. The fit is very good at low frequency. The two curves cross just beyond the 
highest frequency shown, and large deviations occur there. 

Returning to Fig. 7, the short-dashed-long-dashed curve is given for T/Tc = 10 (290 
K). The width of the imaginary part of the conductivity has increased drastically over its 
value for T/T^ = 1. This increase is due entirely to the increased inelastic scattering at 
elevated temperatures. The total width, of course, is a combination of the inelastic and 
elastic contribution; this is clear from a comparison of the T/T^ = 1 curves of Figs. 5 and 
7. In Fig. 5, the width of the imaginary part of the conductivity is very narrow because 
it is only due to inelastic scattering while in Fig. 7 it is almost entirely due to the elastic 
scattering provided by a 1/r = 50 meV impurity term. While one might think that a 
lower bound for this width is 50 meV, one must recall that the effective elastic scattering 



rate is renormalized by the electron-phonon interaction. It is well known ||30|j3^ that the 
electron-phonon interaction renormalizes 1/r by a factor 1/(1 + A), which is about 0.5 in 
our calculations. This is why the effective width in Fig. 7 appears to be smaller than 50 
meV, and close to 25 meV. 



V. THE SPECIFIC CASE OF Bai^^K^BiO^ 



Puchkov et al. p3|j3^ have obtained experimental results for the conductivity in samples 



of Bai^xKxBiO^ with varying x in the superconducting region of the phase diagram. Based 



on an analysis of the real part of the conductivity, the authors of Ref. BS] conclude, contrary 



to the conclusions of many other workers P3-|2H] , that this material cannot be a conventional 
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electron-phonon superconductor with coupling constant value A ~ 1. This conclusion was 



reinforced by the more precise microscopic analysis provided by two of us [|T^. Here, we 
extend our theoretical analysis to the imaginary part. In Fig. 9 we plot the quantity z/cr2(z/) 
as a function of frequency at two different temperatures, one in the superconducting state, 
and one in the normal state. The experimental spectra were obtained using a Kramers- 
Kronig analysis of the reflectivity spectra measured from 25 cm~^ to 40000 cm^l as was 
reported elsewhere |^3|. To improve the accuracy of the analysis, ellipsometric measurements 
were performed on the same crystal at frequencies from 10000 to 50000 cm~^, and the high 
frequency reflectivity approximation was chosen in such a way that the calculated high 
frequency optical conductivity was in agreement with the direct ellipsometric results. Note 
that we have used the data at 1500 cm~^ (185 meV) to normalize the curves. 

The most striking feature of the data in Fig. 9 is the well-defined minimum which occurs 
at approximately 12 meV. In light of the discussion thus far we view this result as unequivocal 
evidence of an s-wave order parameter in the Bai^^K^BiO^ system. Similar results for 
other dopant concentrations (and hence T^'s) are also obtained [Q. This conclusion is 
reinforced by complementary work by Jiang and Carbotte wherein they investigate 
the behaviour of the imaginary part of the conductivity for a superconductor with order 
parameter of d-wave symmetry. They find that no minimum occurs in va2{i') at 2 A, in 
agreement with experimental results on YBaCu^Oj^x, which is generally believed to be a 
d-wave superconductor. This is what is expected since the real part of the conductivity does 
not have an abrupt onset at 2A, which we have argued accompanies the minimum in the 
imaginary part via a Kramers-Kronig transform. 

Given the normalization of the experimental data discussed above. Fig. 2 makes it 
clear that a considerable amount of impurity scattering is required to achieve a qualitative 
agreement between theory and experiment. This observation is also consistent with the 
analysis of the real part of the conductivity, although in that instance a "mid-infrared 
contribution" of unknown origin had to be subtracted from the data first. Then Fig. (6b) 
shows that for A of order unity, the large amount of inelastic scattering present at relatively 

13 



high frequencies (~ 25 — 50 meV) will reduce the imaginary part of the conductivity, in 
disagreement with experiment. This disagreement will be even more significant at higher 
temperatures. 

Theoretical curves are also plotted in Fig. 9 to more clearly illustrate these remarks. As 
a "good" representative fit we show results at the two temperatures using a reduced a^Fi^uj) 
for Bai^xKxBiOi with A = 0.2 (dashed curves). We have used 1/r = 25 meV, a value which 



is consistent with that obtained experimentally from ||33| , M . It is clear that the fit 

is quantitatively good. For comparison we also show the result using the full spectrum, i.e. 
with A = 1. It is obvious that such a fit is poor, and given our previous results, the reader 
can appreciate that no amount of parameter adjustment will yield good agreement with the 
data at both temperatures while retaining A of order unity. This general result is, of course, 
consistent with the conclusion inferred from the real part of the conductivity |33|JT3[1 . The 
analysis here, however, has the advantage that no data subtraction is required before the 
analysis. It has the disadvantage that a high frequency scale is required to normalize the 
experimental data shown in Fig. 9. 
VI. CONCLUSIONS 

We have studied the frequency dependence of the imaginary part of the conductivity 
i/(j2(z/) for an s-wave superconductor. Universal BCS results are first established which 
show a sharp cusp-like minimum in va2{i') at twice the gap value. This provides a very 
good method for determining the gap value and its temperature dependence. As the rate of 
elastic impurity scattering is increased, the gap structure becomes progressively smeared and 
shifts to higher frequency. This shift is completely uncorrelated with the gap value since the 
gap is independent of elastic impurity scattering in an isotropic s-wave superconductor. In 
the presence of impurity scattering strong coupling effects modify the results by smearing the 
gap minimum, and renormalizing the low frequency behaviour through a renormalization 
1/f ^ IIMI- addition, structure is apparent in the phonon region. These are all 



quantitative changes. In contrast, qualitative changes occur in the clean limit, since then 
strong coupling effects become the only source for electron scattering. In this case, there is 
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no structure in 1^(72(1^) at twice the gap value because absorption across the gap can only 
proceed through a phonon- assisted mechanism needed to absorb the necessary momentum 
and energy. This means that the real part of the conductivity will become non-zero at 
a frequency of twice the gap plus the lowest phonon energy in the system. The phonon 
spectrum generally extends to zero frequency but the density of states for the phonons 
usually has an lo^ dependence at low frequency. This leads to a rather gradual increase 
in absorption and no sharp structure in va2{y)- When impurities are added, a qualitative 
change occurs because impurity-assisted absorption can now set in quite abruptly at 2 A. 
This manifests itself in the imaginary part of the conductivity as a sharp cusp-like minimum 
in vGiiy). 

Application of our calculations to the specific case of Bax^^K^BiO-i, allows us to conclude, 
from a consideration of data on the imaginary part of the conductivity alone, that this 
system is an s-wave superconductor. On the other hand further analysis suggests that the 
superconductivity is not phonon-mediated as others have concluded. The data is consistent 
with an electron- phonon mass enhancement parameter of about A ~ 0.2. Larger values of 
order unity suggested in the literature lead to much too large a scattering rate at room 
temperature. 

In summary we believe that measurement of the imaginary part of the conductivity is a 
powerful alternative to determine both the symmetry and the value of the superconducting 
order parameter. Used with normal state measurements, such a probe also allows us to infer 
the strength of the inelastic scattering present in a particular superconductor. In contrast, 
in a d-wave superconductor, there will be no gap minimum at 2A because in such a case 
there is an entire distribution of gaps, starting at value zero. The conductivity averages over 
these gaps, and this leads to smooth behaviour in va2{y)- 
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Figure Captions 

• Fig. 1. The quantity 1^(72(1^) (imaginary part of the conductivity times the frequency 
in units of ne^/m) vs. frequency for various impurity scattering rates, (a) l/r — 1, (b) 
1/r = 5, and (c) l/r = 25 meV. These results are calculated in the BCS limit with 
Tc — 29 K. The zero temperature gap value is given by = 3.53 (Aq = 8.8 meV). 
The various curves correspond to different temperatures: T/Tc = 0.1 (long-dashed), 
T /Tc = 0.5 (long-dashed-short-dashed), T/Tc = 0.7 (long-dashed-dotted), T/Tc = 0.8 
(short-dashed), T/T^ = 0.9 (dotted), and T/T^ = 0.995 (solid). 

• Fig. 2. i'a2{i') vs. frequency for various impurity scattering rates. These results are in 
the BCS hmit with T/T^ = 0.1. The scattering rates are l/r = 1 meV (sohd), l/r = 5 
meV (dotted), l/r = 10 meV (dashed), l/r = 25 meV (dot-long-dashed), l/r = 50 
meV (short-dashed-long-dashed), and l/r = 100 meV (long-dashed). 

• Fig. 3. i>a2{iy) vs reduced temperature T/Tc in the BCS limit, plotted for various 
frequencies. The results are for (a) l/r = 1 and (b) l/r = 25 meV. The curves 
correspond to the frequencies u — 0.2 meV (solid), v — 1.0 meV (dotted), v — 2.5 
meV (dashed), z/ = 5.0 meV (dot-dashed), and u = 10 meV (short-dashed-long- 
dashed). The effect of finite frequency is more significant in the clean limit. 

• Fig. 4. 1/(72(1^) vs the impurity scattering rate l/r, for various frequencies. The results 
are in the BCS limit, with T/Tc — 0.1. The curves correspond to u — (sohd), 1 
(dotted), 5 (dashed), and 10 meV (dot-dashed). 

• Fig. 5. 1^172(1^) vs. frequency as calculated using an electron-phonon spectrum with 
A = 1 and = 29 K as described in the text. In (a) we use l/r = (clean hmit) 
and plot the results for various temperatures, as indicated. The solid (dotted) curve 
is for the superconducting (normal) state at T/Tc — 0.99. Clearly there is very little 
difference throughout the frequency scale shown. Note the lack of a cusp-Hke minimum 
ai u — 2Ao ~ 10.4 meV. In (b) we focus on the low temperature results {T/Tc — 0.2) 
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and show the effect of a small amount of impurity scattering. The cusp-like minimum 
at 2 A is clearly present. 

Fig. 6. 

i/(72(i^) VS. frequency for various coupling strengths, as indicated. The results are for 
low temperature T/Tc — 0.3 and (a) relatively clean, 1/t — 2 meV, and (b) 1/r = 25 
meV. The Coulomb repulsion yU* has been adjusted to keep Tc fixed while the coupling 
strength has changed (simply by scaling the a^F{uj) spectrum). Increased coupling 
strength suppresses v(T2{y) and broadens the minimum at 2A. Note that 2A increases 
slightly as the coupling strength is increased. 

Fig. 7. z/(J2(z/) vs. frequency for 1/r = 50 meV. The other parameters are the same 
as used in Fig. 5. The various curves are for difi^erent temperatures, T/T^ = 0.99 
(solid — superconducting, dotted — normal), T /Tc — 0.5 (dashed), and T /Tc — 0.2 
(long-dashed). Both the first pair and the second pair of curves differ from one another 
only at very low frequency. Also shown is the normal state result at T/Tg = 10. 

Fig. 8. 1/(72 (z/) vs. frequency for l/r = (clean limit) in the normal state, using the 
same parameters for the spectrum as in Fig. 5. The solid curves are results of complete 
Eliashberg calculations whereas the dotted curves are for the derived low frequency 
Drude fits to the complete results. Curves are labelled by their reduced temperatures. 

Fig. 9. The measured va2{i') vs. frequency at T = 9 K and at T = 300 K (solid 
curves). Also shown are the theoretical fits, using the a^F{uj) described in the text, 
but with A = 0.2 (dashed curves). Tc is kept fixed to the experimental value with a 
negative /x*. Finally, theoretical fits are also shown with A = 1 (dotted curves). The 
latter curves are clearly incompatible with the experimental results. 
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